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We have made a new evaluation of the A decay width in nuclear matter 
within the Propagator Method. Through the Local Density Approximation 
^NJ . it is possible to obtain results in finite nuclei. We have also studied the 

dependence of the widths on the NN and AN short range correlations. Using 
reasonable values for the parameters that control these correlations, as well 
as realistic nuclear densities and A wave functions, we reproduce, for the first 
time, the experimental non-mesonic widths in a wide range of mass numbers 
(from medium to heavy hypernuclei) . 
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I. INTRODUCTION 

A hypernucleus is a bound system made of by neutrons, protons and one or more hyper- 
ons. Among these strange nuclei, those which contain one A hyperon are the most stable 
with respect to the strong interaction and they are the subject of this paper. The study of 
hypernuclear physics may help in understanding some present problems related, for instance, 
to some aspects of weak interactions in nuclei, or to the origin of the spin-orbit interaction 
in nuclei. Besides, it is a good instrument to study the role of quark degrees of freedom 
in the hadron-hadron interactions at short distances and the renormalization properties of 
pions in the nuclear medium. 

Nowadays we know some important features of the YN interaction fl|]. For example, at 
intermediate distances the strong AN interaction is weaker than the NN one, and its spin- 
orbit term is very small. Moreover, the former has a smaller range than the NN one. From 
the study on mesonic decays of light hypernuclei we have evidence for strongly repulsive 
cores in the A A" interaction at short distances 0, which automatically appears in quark 
based models P,|J. These characteristics of the AA" interaction are important, as we will 
see, for the evaluation of the decay rates of A-hypernuclei. 

The most interesting hypernuclear decays are those involving weak processes, which 
directly concern the hyperon. The weak decay of hypernuclei occurs via two channels: the so 
called mesonic channel (A — > irN) and the non-mesonic one, in which the pion emitted from 
the weak hadronic vertex is absorbed by one or more nucleons in the medium (AN — > NN, 
ANN — > NNN, etc. ). Obviously, the non-mesonic processes can also be mediated by the 
exchange of more massive mesons than the pion. The non-mesonic decay is only possible in 
nuclei and, nowadays, the study of the hypernuclear decay is the only practical way to get 
information on the weak process AN — > NN, especially on its parity conserving part. In 
fact, there are not experimental observations for this interaction using lambda beams. It is, 
however, under study the inverse reaction pn — > Ap at COSY || and RCNP |J. 

The free A decay is compatible with the AI = 1/2 isospin rule, which is also valid for 
the decay of other hyperons and for kaons (namely in non-leptonic strangeness changing 
processes). This rule is based on the experimental observation that the A — > ir~p decay 
rate is twice the A — > n°n one, but it is not yet understood on theoretical grounds. From 
theoretical calculations like the one in ref. @ and from experimental measurements |§ there 
is some evidence that the AI = 1/2 rule is broken in nuclear mesonic decay. However, this 
is essentially due to shell effects and might not be directly related to the weak process. A 
recent estimate of AI = 3/2 contributions to the AN — > NN reaction || found moderate 
effects on the hypernuclear decay rates. In the present calculation of the decay rates in 
nuclei we will assume this rule as valid. The momentum of the final nucleon in A — > nN is 
about 100 MeV for A at rest, so this process is suppressed by the Pauli principle in nuclei 
(particularly in heavy systems). It is strictly forbidden in infinite nuclear matter (where 
kp ~ 270 MeV), but in finite nuclei it can occur because of three important effects: 1) in 
nuclei the hyperon has a momentum distribution that allows larger momenta for the final 



nucleoli, 2) the final pion feels an attraction by the medium such that for fixed momentum it 
has a smaller energy than the free one and consequently, due to energy conservation, the final 
nucleon again has more chance to come out above the Fermi surface, and 3) on the nuclear 
surface the local Fermi momentum is smaller than k° F and favours the decay. Nevertheless, 
the mesonic width decreases fastly as the mass number A of the hypernucleus increases ]7| . 
From the study of the mesonic channel it could be possible to extract important information 
on the pion-nucleus optical potential, which we do not know today in a complete form. In 
fact, the mesonic rate is very sensitive to the pion self-energy in the medium M. 

The final nucleons in the non-mesonic process AN — > NN emerge with large momenta 
(~ 420 MeV), so this decay is not forbidden by the Pauli principle. On the contrary, apart 
from very light hypernuclei (the s-shell ones) , it dominates over the mesonic decay. The non- 
mesonic channel is characterized by large momentum transfers, so that the details of the 
nuclear structure do not have a substantial influence while the NN and AiV Short Range 
Correlations (SRC) turn out to be very important. There is an anticorrelation between 
mesonic and non-mesonic decay modes such that the total lifetime is quite stable from light 
to heavy hypernuclei [§,0]: r exp = (0.5 -f- 1) r/ ree - 

Nowadays, the main problem concerning the weak decay rates is to reproduce the ex- 
perimental values for the ratio T n /T p between the neutron and the proton induced widths 
An — > nn and Ap — > np. The theoretical calculations underestimate the experimental data 
for all the considered hypernuclei PJ9|JTlHT4]| : 




Exp / -p \ Exp 

0.5 < ^ <2. (1.1) 



In the One Pion Exchange (OPE) approximation the values for this ratio are 0.1 -f- 0.2. 
On the other hand the OPE model has been able to reproduce the 1-body stimulated non- 



mesonic rates Ti = T n + T p for light and medium hypernuclei f|n|,|lj,|lj] . In order to solve 
this problem many attempts have been made up to now, but without success. Among these 
we recall the inclusion in the AA" — > NN transition potential of mesons heavier than the 
pion |TT|JT3| , p!4]| , the inclusion of interaction terms that violate the AI = 1/2 rule §| and 
the description of the short range baryon-baryon interaction in terms of quark degrees of 
freedom |12| . This last calculation is the only one that has found a consistent (but not 
sufficient) increase of the neutron to proton ratio with respect to the OPE one. However, 
this calculation is only made for s-shell hypernuclei and their effective quark-lagrangian 
does not reproduce the experimental ratio between the AI = 1/2 and AI = 3/2 transition 
amplitudes for the A free decay. 

The analysis of the ratio r n /r p is influenced by the 2-nucleon induced process ANN — > 
NNN. By assuming that the meson produced in the weak vertex is mainly absorbed by 
a neutron-proton strongly correlated pair, the 3-body process turns out to be Anp —> nnp, 
so that a considerable fraction of the measured neutrons could come from this channel and 
not only from the An — > nn and Ap — > np ones. In this way it might be possible to explain 
the large experimental r n /r p ratios, which originally have been analyzed without taking 



into account the 2-body stimulated process. Nevertheless, the situation is far from being 
clear and simple. The new non-mesonic mode was introduced in ref. |15] and its calculation 
was improved in ref. [16], where the authors found that the inclusion of the new channel 



would lead to larger values of the T n /T p ratios extracted from the experiment, somehow 
more in disagreement with theoretical estimates. However, in the hypothesis that only two 
nucleons from the 3-body decay are detected, the reanalysis of the experimental data would 
lead to smaller ratios ||. These observations show that T n /T p is sensitive to the energy 
spectra of the emitted nucleons, whose calculation also requires a careful treatment of the 
nucleon Final State Interaction. In ref. [|r7| the energy distributions were calculated using 
a Monte Carlo simulation to describe the final state interactions. A direct comparison of 
those spectra with the experimental ones favours T n /T p values around 2 -=- 3 (or higher), in 
disagreement with the OPE predictions. However, it was also pointed out the convenience of 
measuring the number of protons per decay event. This observable, which can be measured 
from delayed fission events in the decay of heavy hypernuclei, gives a more reliable ratio 
T n /T p and is less sensitive to details of the Monte Carlo simulation determining the final 
shape of the spectra. 

In this paper we present a new evaluation of the decay rates for medium to heavy 
hypernuclei based on the Propagator Method of ref. |18j, which allows a unified treatment 



of all the decay channels. The parameters of the model are adjusted to reproduce the non- 
mesonic width of \ 2 C and the decay rates of heavier hypernuclei are predicted. We also 
discuss how the new model affects the energy spectrum of the emitted nucleons, in the hope 
of obtaining a ratio T n /T p more in agreement with the experimental observation. 

The paper is organized as follows. In Sec. [TI|we present the model used for the calculation 
of the decay rates. Our results are presented and discussed in Sec. [ITT[ We first study 
the sensitivity of the decay rates to the parameters defining the NN and AN short range 
correlations as well as to the nuclear density and A wave functions. We then obtain the 
widths for various hypernuclei and discuss the energy distribution of the nucleons from the 



weak decays. Our conclusions are given in Sec. [TV 



II. PROPAGATOR METHOD 



The A decay in nuclear systems can be studied in Random Phase Approximation (RPA) 
using the Propagator Method |TB|,[rS|. This technique provides a unified picture of the dif- 
ferent decay channels and it is equivalent to the standard Wave Function Method (WFM) 
|n§ , used by other authors in refs. @,P,|n],|l3|,|l4[] . The calculation of the widths is usually 



performed in nuclear matter, and then extended to finite nuclei via the Local Density Ap- 
proximation (LDA). For the calculation of the mesonic rates the WFM is more reliable than 
the propagator method in LDA since this channel is rather sensitive to the shell structure of 
the hypernucleus, given the small energies involved. Moreover, it is advisable to avoid the 
use of the LDA for very light systems and we will make the calculation starting from \ 2 C. 
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FIG. 1. A self energy in nuclear matter 



On the other hand, the propagator method in LDA offers the possibility of calculations over 
a broad range of mass numbers. 



The method was introduced in ref. [|jj| and we briefly summarize it here for clarity. The 
A — > it N effective lagrangian is: 



C 



AnN 



Gmlip N (A + 5 75 )r ■ 0> A + h.c. 



(2.1) 



where the values of the weak coupling constants G ~ 2.211 • 10~ 7 /m 2 , A = 1.06, B = —7.10 
are fixed on the free A decay. The constants A and B determine the strengths of the parity 
violating and parity conserving A — > itN amplitudes respectively. In order to enforce the 
A/ = 1/2 rule, in eq. ( j2.1|) the hyperon is assumed to be an isospin spurion with I z = —1/2. 
To calculate the A width in nuclear matter we start with the imaginary part of the A 
self-energy: 



-2Im£ A . 



(2.2) 



By the use of Feynman rules, from fig. [I] it is easy to obtain the A self-energy in the following 
form: 



£ A (A-) = -ii(Gm*) 2 l-^-\s 2 + ^q 2 )F*{q)G N {k-q)GM)- 



(2rr)< 



mi 



(2.3) 



Here, S = A, P = m^B /2rriN, while the nucleon and pion propagators in nuclear matter 
are respectively: 



Gn(p) 



0(1 P I -kp) 



+ 



9(k F - | p |) 



p -E N (p)-V N + ie p -E N (p) - V N - ie 



(2.4) 



and: 

G -(<?) = 2 \ w v ( 2 - 5 ) 

q$-q 2 -ml - S* (g) 

In the above, p = (po,P) an d 9 — (%,Q) denote four-vectors, kp is the Fermi momentum, 
E N is the nucleon total free energy, Vjq is the nucleon binding energy, and S* is the pion 
proper self-energy in nuclear matter. Moreover, in eq. (|2.3|) we have included a monopole 
form factor for the 7rAiV vertex: 

A 2 _ jy.'i 

(the same is used for the nNN strong vertex), with cut-off A„- = 1.2 GeV. In fig. ^| we 
show the lowest order Feynman graphs for the A self-energy in nuclear matter. Diagram (a) 
represents the bare self-energy term, including the effects of Pauli principle and of binding 
on the intermediate nucleon. In (b) and (c) the pion couples to a particle-hole (p-h) and 
a A-h pairs, respectively. Diagram (d) is an insertion of s-wave pion self-energy at lowest 
order. In diagram (e) we show a 2p-2h excitation coupled to the pion through s-wave tcN 
interactions. Other 2p-2h excitations, coupled in p-wave, are shown in (f), while (g) is a 
RPA iteration of diagram (b). It is possible to evaluate the integral over q in (|2.3| ), and the 
A self-energy (eq. ( [2.2|) ) in nuclear matter becomes ||18|| : 

dq 
xlma;(g) \ qo=ko _E N (k- q )-v N , (2.7) 



T A (fc, p) = -6(Gml) 2 J -7^0(\ k-q\ -k F )9(k - E N (k - q) - V N ) 



where 



a{q) 



S 2 (q)U(q) 







V L (q)U(q) l-V T (q)U(q) 



(2.8) 



In eq. (|2.7|) the first 9 function forbids intermediate nucleon momenta (see fig. |I]) smaller 
than the Fermi momentum and the second one requires the pion energy q to be positive. 
Moreover, the A energy, k$ = E\(k) + V\, contains a binding term. The pion lines of fig. |2| 
have been replaced in eq. fl2.8p by the interactions S, Pl, Pt,Vl, Vt, which include n and p 
exchange modulated by the effect of short range correlations and whose espressions are given 
in the Appendix. The functions Vl and Vt represent the (strong) p-h interaction, including 
a short range Landau parameter g', while S, Pl and Pt correspond to the lines connecting 
weak and strong hadronic vertices and contain another short range Landau parameter g' A . 
Furthermore, in eq. ( |2.8|) : 

Gl(q) = \ -, (2.9) 

q 2 -q 2 - m 2 n 
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FIG. 2. Lowest order terms for the A self-energy in nuclear matter. The meaning of the various 
diagramms is explained in the text. 



is the free pion propagator, while U(q) contains the Lindhard functions for p-h and A-h 



excitations ]20] and also accounts for 2p-2h excitations: 

U(q) = U ph {q) + U Ah (q) + U 2p2h {q). (2.10) 



It appears in eq. (|2.S| ) within the standard RPA expression. Eq. ( |2.7|) depends explicitly 
and through U(q) on the nuclear matter density p = 2k F /3ir 2 . The Lindhard functions U p h, 
U/\h are normalized as in ref. pi] and U2 P 2h is evaluated as in |j|, that is, by calculating 
the available phase space for 2p-2h excitations and by taking into account the experimental 
data on pionic atoms. U(q) is related to the pion proper self-energy through: 



^q 2 Fl{q)U{q) 



K(i) = ^*(q) + 4 S) *(?), 4 P) *0?) = -^ > (2.113 

^9L(q)U(q) 



mi 



where the Landau function <7l(<?) is given in the Appendix [see eq. (|A12|) 1, and E^* is 



the s-wave part of the self-energy We will use the parametrization of ref. p2| : S^*- ) *(g) = 
— 47r(l + m 7r /mAr)6 p, with 60 = — 0.0285/m„-. The function Sj 8 '* is real (constant and 
positive), therefore it contributes only to the mesonic decay (diagram (d) in fig. |2| is the 
relative lowest order). On the contrary, the p-wave self-energy E^ * is complex and attractive 
(that is, ReE<f)*(g) < 0). 

The decay widths in finite nuclei are obtained in LDA. In this approximation, the Fermi 
momentum becomes r-dependent (that is, a local Fermi sea of nucleons is introduced) and 
related again to the nuclear density by: 

k F (r) = [-^p(r)j . (2.12) 

Besides, the nucleon binding potential Vjv also becomes r-dependent in LDA. In Thomas- 
Fermi approximation we assume: 

e F {r) + V N {r) = ^- + V N (r) = 0. (2.13) 

ZrriN 



For the A binding energy we use instead the experimental value p3| , |24|| . With these pre 



scriptions we can then evaluate the decay width in finite nuclei through the relation: 

r A (fc) = J dr\ Mr) I 2 r A [k,p(r)], (2.14) 



where ip\ is the A wave function and I\[fc, p(r)] is given by eqs. (p.7|) , (|2.8|) . This decay 



rate is valid for fixed A momentum k. A further average over the A momentum distribution 
gives the total width: 

r dk\Mk)\ 2 T A (k), (2.15) 



which can be compared with experimental results. 

The propagator method provides a unified picture of the decay widths. The imaginary 
part of a self-energy diagram requires placing simultaneously on-shell the particles of the 
considered intermediate state. For instance, diagram (b) in fig. |^ has two sources of imagi- 
nary part. One comes from cut 1, where the nucleon and the pion are placed on-shell. This 
term contributes to the mesonic channel: the final pion interacts with the medium through 
a p-h excitation and then escapes from the nucleus. Diagram (b) and further iterations lead 
to a renormalization of the pion in the medium which increases the mesonic rate by about 
two orders of magnitude in heavy nuclei JT8|. The cut 2 in fig. 0(b) place a nucleon and a 



p-h pair on shell, so it is the lowest order contribution to the physical process AN — > NN. 
The mesonic width r^ is calculated from 

a(q) = ls 2 + ^q 2 \G n (q), (2.16) 

by omitting ImE* in G w , namely by replacing: 

lmG n (q) - -n5(q 2 - q 2 - m* - ReEfa))- (2.17) 

The one-body induced non- mesonic decay rate T 1 is obtained by substituting in eqs. fl2.7Q, 



l-V LiT (q)U(q) \l-V L>T (q)U(q)r l ' ' 

that is by omitting the imaginary parts of U&h and U2 P 2h in the numerator. Indeed Ihi[/a/i 
accounts for the A — > itN decay width, thus representing a contribution to the mesonic decay. 
There is no overlap between ]mU P h(q) and the pole go = u(o) m eq. ( |2.17| ), so the separation 



of the mesonic and 2-body non-mesonic channels is unambiguous. The renormalized pion 
pole in eq. fl2.16| ) is given by the dispersion relation: 



u 2 (q) -q 2 -ml- ReZ;[uj(q), q] = 0, (2.19) 

with the constraint: 

q = k -E N (k-q)-V N . (2.20) 

At the pion pole ImU2 P 2h 7^ 0, thus the 2-body induced non-mesonic width 1^ cannot be 
calculated using the prescription ( |2.18|) with U2 P 2h instead of U p h in the numerator of the 



r. h. s. Part of the decay rate calculated in this way it is due to excitations of the renormalized 
pion and contributes to Tm- The 3-body non-mesonic rate is then calculated by subtracting 
T M and I\ from the total rate T TO t, which we get via the full espression for a [Eq. (|2.8|) ]. 



III. RESULTS AND DISCUSSION 

Let us now discuss the numerical results one can obtain from the above illustrated formal- 
ism. We shall first study the influence of short range correlations and the A wave function on 
the decay width of \ 2 C, which will be used as a testing ground for the theoretical framework 
in order to fix the parameters of our model. We will then obtain the decay widths of heavier 
hypernuclei and we will explore whether the refined model influences the energy distribution 
of the emitted particles, following the Monte Carlo procedure of ref. [T7\ . 



In order to evaluate the width ( [2.14| ) in LDA one needs to specify the nuclear density 



and the wave function for the A. The former is assumed to be a Fermi distribution: 

.4 



i ^ P^ A ) 



1 _|_ e [r-R(A)]/a 



po(A) i^(A){i+[^n 



(3.1) 



with radius R(A) = 1.12A 1/3 - 0.86A _1/3 fm and thickness a = 0.52 fm. The A wave 
function is obtained from a Wood-Saxon (W-S) well which exactly reproduces the first two 
single particle eigenvalues (s and p levels) measured in A-hypernuclei. 

A. Short range correlations and A wave function 

A crucial ingredient in the calculation of the decay widths is the short range part of the 
AW and AA" interactions. They are expressed by the functions gL,T{l) and g^ T (q), which are 
reported in the Appendix and contain the Landau parameters g' and g' A . No experimental 
information is available on g' A , while many constraints have been set on g', for example by 
the well known quenching of the Gamow- Teller resonance. Realistic values of g', within the 
framework of the ring approximation, are in the range 0.6 -r 0.7 p5|. However, in the present 



context g' correlates not only p-h pairs but also p-h with 2p-2h states. In order to fix these 
parameters we shall compare our calculations with the experimental non-mesonic width of 

12 r\ 

A ^- 

In fig. |H we see how the total non-mesonic width for carbon depends on the Landau 
parameters. The rate decreases as g' increases. This characteristic is well established in RPA. 
Moreover, fixing g', there is a minimum for g' A ~ 0.4 (almost independent of the value of g'). 
This is due to the fact that for g A -C 0.4 the longitudinal p-wave contribution in eq. ( |2.8| ) 
dominates over the transverse one and the opposite occurs for g' A ^> 0.4. We also remind 
that the s-wave interactions are independent of g A [eq. flAlC )]. Moreover, the longitudinal 



p-wave AA" — > NN interaction [eq. ( |A"8|) ] contains the pion exchange plus SRC, while the 
transverse p-wave AA" — > AW interaction [eq. (|A9|) 1 only contains repulsive correlations, so 
with increasing g' A the p-wave longitudinal contribution to the width decreases, while the p- 
wave transverse part increases. From fig. |3] we see that there is a broad range of choices of g' 
and g' A values which fit the experimental band. The latter represents the non-mesonic decay 
width which is compatible with both the BNL |2(J and KEK |2TJ experiments. One should 

10 




g'A 

FIG. 3. Dependence of the non-mesonic width on the Landau parameters g' and g'^ for ^ 2 C. 
The experimental value from BNL |2(| (KEK |[27l ) lies in between the horizontal solid (dashed) 
lines. 
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TABLE I. Wave function sensitivity 



12 n 


H.O. 


Dover 


New W-S 


Microscopic 


BNL JH 


KEK [27| 


KEK New [3^,31] 


^M 


0.26 


0.25 


0.25 


0.25 


0.11 ±0.27 


0.36 ±0.15 


>0.11 


Ti 


0.78 


0.77 


0.82 


0.69 








r 2 


0.15 


0.15 


0.16 


0.13 








^NM 


0.93 


0.92 


0.98 


0.81 


1.14 ±0.20 


0.89 ±0.18 


< 1.03 


Ttot 


1.19 


1.17 


1.23 


1.06 


1.25 ±0.18 


1.25 ±0.18 


1.14 ±0.08 



notice that the theoretical curves reported in fig. |3] contain the contribution of the 3-body 
process; should the latter be neglected (ring approximation) then one could get equivalent 
results with g' values smaller than the ones reported in the figure (tipically Ag' ~ —0.1). 
The phenomenology of the (e, e') quasi-elastic scattering suggests, in ring approximation, g' 
values of the order of 0.7. Here, by taking into account also 2p — 2h contributions, we shall 
use the "equivalent" value g' = 0.8, together with g' A = 0.4. 

We note that the values used in ref. |16[], namely g' = 0.615 and g' A = 0.62, would yield 
r"i = 1.26 and T 2 = 0.25, adding to a non-mesonic width T^m — 1-51, which is 50% larger 
than the experimental one. Thus the analysis performed here shows that the present data 
for ^ 2 C favour a somewhat different but still reasonable g' value. 

We shall now illustrate the sensitivity of our calculation to the A wave function in \ 2 C In 
addition to the W-S that reproduces the s and p levels, we also use a harmonic oscillator wave 



function with an "empirical" frequency uo [[2^,[24|], again obtained from the s — p energy shift, 
the W-S wave function of Dover et al. |28| and the microscopic wave function calculated 

29| . The results are 



from a non-local self-energy using a realistic YN interaction in ref. 

shown in table |, where they are compared with the experimental data from BNL |26j and 

KEK PSil. 



By construction, the chosen g' and g' A reproduce the experimental decay widths using 
the W-S wave function which gives the right s and p levels in ^ 2 C We note that it is possible 



to generate the microscopic wave function of ref. |29[] for carbon via a local hyperon-nucleus 
W-S potential with radius 2.92 fm and depth —23 MeV. Although this potential reproduces 
fairly well the experimental s-level for the A in \ 2 C, it does not reproduce the p-level. In this 
work we prefer to use a completely phenomenological A-nucleus potential that can easily 
be extended to heavier nuclei and reproduces the experimental A single particle levels as 
well as possible. Except for s-shell hypernuclei, where experimental data require A-nucleus 
potentials with a repulsive core at short distances 0, the A binding energies have been well 
reproduced by W-S potentials. We thus use a W-S potential with fixed diffuseness (a = 0.6 
fm) and adjust the radius and depth to reproduce the s and p A-levels. The paramaters of 
the potential for carbon are R = 2.27 fm and Vq = —32 MeV. 

To analyse the results of table |, we note that the microscopic wave function is substan- 
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TABLE II. W-S parameters 



A Z 


i? (fm) 


Vq (MeV) 


A u 


2.27 


-32.0 


fSi 


3.33 


-29.5 


A Ca 


4.07 


-28.0 


^ 6 Fe 


4.21 


-29.0 


89v 
A Y 


5.07 


-28.5 


A 39 La 


6.81 


-27.5 


208p b 


5.65 


-32.0 



tially more extended than all the other wave functions used in the present study. The Dover 
parameters 28|| , namely R = 2.71 fm, Vq = —28 MeV, give rise to a A wave function that is 



somewhat more extended than the new W-S one but is very similar to that obtained from 
a harmonic oscillator with a frequency of 10.9 MeV, adjusted to the s — p energy shift in 
carbon. Consequently, the non-mesonic width from the Dover's wave function is very similar 
to the one obtained from the harmonic oscillator and slightly smaller than the new W-S one. 
The microscopic wave-function predicts the smallest non-mesonic widths due to the more 
extended A wave-function, which explores regions of lower density and thus has a smaller 
probability of interacting with one or more nucleons. From table | we also see that, against 
intuition, the mesonic width is quite insensitive to the A wave function. On this point we 
remind that for fixed A-momentum, the more extended is the wave function in r-space, 
the larger is the mesonic width, since the Pauli blocking effects on the emitted nucleon are 
reduced. But, when we make the integral over the A- momenta in LDA (eq. (|2.15 )), to more 



extended wave functions in r-space correspond less extended momentum distributions which 
tend to decrease the mesonic width. The two effects tend to cancel each other and Tm is 
insensitive to the different wave functions used in the calculation. In summary, different 
(but realistic) A wave functions give rise to total decay widths which may differ at most by 

15%. 

B. Decay widths of medium-heavy hypernuclei 

Using the new W-S wave functions and the Landau parameters g' = 0.8, g'^ = 0.4 we 
have extended the calculation to heavier hypernuclei. We note that, in order to reproduce 
the experimental s and p levels for the hyperon we must use potentials with nearly constant 
depth, around 28 -f- 32 MeV, from medium to heavy hypernuclei (radii and depths of the 
used W-S potentials are quoted in table [TI|). 



Our results are shown in table III. We observe that the mesonic rate rapidly vanishes 



by increasing the mass number A. This is well known and it is related to the decreasing 



13 



TABLE III. Decay rates 



A Z 


Ta/ 


Ti 


r 2 


Ttot 


A u 


0.25 


0.82 


0.16 


1.23 


fSi 


0.07 


1.02 


0.21 


1.30 


A Ca 


0.03 


1.05 


0.21 


1.29 


^ 6 Fe 


0.01 


1.12 


0.21 


1.35 


89v 
A Y 


6 • 10~ 3 


1.16 


0.22 


1.38 


139t „ 

A 


6 • 10~ 3 


1.14 


0.18 


1.33 


208p b 


1 • io- 4 


1.21 


0.19 


1.40 



phase space allowed for the mesonic channel, and to smaller overlaps between the A wave 
function ip\ and the nuclear surface, as A increases. The 2-body induced decay is rather 
independent of the hypernuclear dimension and it is about 15% of the total width. Previous 
works fll~5| , |16|| gave more emphasis to this new channel, without, however, reproducing the 
experimental results. The total width is also nearly constant with A, as we already know 
from the experiment. In fig. ^ we compare the results from table |TT| with recent (after 1990) 
experimental data for non-mesonic decay |^6, 27, 30, 32-34]. Nevertheless, we remind that the 



data for nuclei from A 8 Si on refer to the total width. However, as can be seen from table [II, 
Fm (a^i) /Tnm( i 8 Si) — 6 • 10 -2 and this ratio rapidly decreases with A. The theoretical 
results are in good agreement with the data (which, on the other hand, have large error bars) 
over the whole hypernuclear mass range explored. Moreover, we also see how the saturation 
of the AA" — > NN interaction in nuclei is well reproduced. 

One of the open problems in the study of weak hypernuclear decays is to understand 
the large experimental value of the ratio T n /T p which most of the present theories fail to 
reproduce. Only the quark model of ref. |12| predicts an enhanced ratio, although it cannot 
describe both mesonic and non-mesonic decays from the same basic quark hamiltonian. 
However, we have to remind that the data for T n /T p have a large uncertainty and they 
have been analyzed without taking into account the 3-body decay mechanism. The study 
of ref. |IIJ showed that, even if the three body reaction is only about 15 % of the total 
decay rate, this mechanism influences the analysis of the data determining the ratio T n /T p . 
The energy spectra of neutrons and protons from the non-mesonic decay mechanisms were 
calculated in ref. [17]. The momentum distributions of the primary nucleons were determined 



from the Propagator Method and a subsequent Monte Carlo simulation was used to account 
for the final state interactions. It was shown that the shape of the proton spectrum was 
sensitive to the ratio r n /r p . In fact, the protons from the three-nucleon mechanism appeared 
mainly at low energies, while those from the two-body process peaked around 75 MeV. Since 
the experimental spectra show a fair amount of protons in the low energy region they would 
favour a relatively larger three-body decay rate or, conversely, a reduced number of protons 
from the two-body process. Consequently, the experimental spectra are compatible with 
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FIG. 4. A decay widths in finite nuclei as a function of the mass number A. 
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FIG. 5. Proton spectrum from the decay of \ 2 C for various values of F n /T p . The experimental 



data are taken from ref. [26 



values for T n /T p around 2 -r 3, in strong contradiction with the present theories. 

The excellent agreement with the experimental decay rates of medium to heavy hyper- 
nuclei obtained here from the Propagator Method with modified parameters, makes it worth 
to explore the predictions for the nucleon spectra. The question is whether this modified 
model affects the momentum distribution of the primary emitted nucleons strongly enough, 
such that good agreement with the experimental proton spectra is obtained without the 
need for very large values for T n /T p . We have thus generated the nucleon spectra from 
the decay of several hypernuclei using the Monte Carlo simulation of ref. [|17| but with our 
modified g', g' A parameters and our more realistic nuclear density and A wave functions. 
The spectra obtained for various values of r n /r p , used as a free parameter in the approach 
of ref. Hi? |, are compared with the BNL experimental data [^6j in fig. [|. We observe that, 
although the non-mesonic widths are smaller by about 35% than those of refs. [|16|,|17| > the 
resulting nucleon spectra, once they are normalized to the same total width, are practically 
identical. The reason is that the ratio r 2 /ri of two-body induced versus one-body induced 
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decay rates is essentially the same in both models (between 0.2 and 0.15 from medium to 
heavy hypernuclei) , and the momentum distributions for the primary emitted protons are 



also very similar. As a consequence, the conclusions drawn in ref. |Ty]] still hold and the 
present calculation would also favour very large values of T n /T p . 

Therefore, the origin of the discrepancy between theory and experiment for the ratio 
T n /T p still needs to be resolved. From the theoretical side, there is still room for improving on 
the numerical simulation of final state interactions. In particular, Coulomb distorsions and 
the evaporating processes need to be incorporated. We think that the evaporating process 
is an important ingredient which increases the nucleon spectra at low energies. Maybe this 
contribution is so important that there is no need for high T n /T p values. On the experimental 



side, although new spectra are now available f27jj3(| , they have not been corrected for energy 



losses inside the target or detector, so a direct comparison with the theoretical predictions 
is not yet possible. Attempts to incorporate these corrections by combining a theoretical 
model for the nucleon rescattering in the nucleus with a simulation of the energy losses in 
the experimental set-up are now being pursued [|35 ]. These efforts call for newer improved 



theoretical models that incorporate those final state interaction effects missing in ref. |T7| . 
On the other hand, a forward step towards a clean extraction of the ratio T n /T p would 
be obtained if the nucleons from the different non-mesonic processes, AiV — > NN and 
ANN — > NNN were disentagled. Through the measurement of the coincident spectra of 
the outgoing nucleons, it could be possible, in the near future, to split the non-mesonic decay 
width into its two components Ti and 1^ [|36| and obtain a cleaner measurement of the ratio 

J- 71/ J- p- 

IV. CONCLUSIONS 

Using the Propagator Method in Local Density Approximation, in this paper we made 
a new evaluation of the A decay widths in nuclei. Special attention has been devoted 
to the study of the NN and AA" short range interactions and realistic nuclear densities 
and A wave functions were used. We have adjusted the parameters that control the short 
range correlations to reproduce the experimental decay widths of \ 2 C. Then, the calculation 
has been extended to heavier hypernuclei, up to 208 Pb. We reproduce for the first time 
the experimental non-mesonic decay widths from medium to heavy A-hypernuclei and the 
saturation of the AA" — > NN interaction is observed. 

The energetic spectra of emitted nucleons calculated using the Propagator Method with 
modified parameters (describing the energy distributions of primary nucleons) and the Monte 
Carlo simulation (accounting for the final state interactions) does not change appreciably 
with respect to those calculated in ref. Hl7fl . The reason is that, in spite of the fact that the 



non-mesonic decay widths Ti and T 2 are sizably reduced (by about 35%) with respect to 
those of ref. |17[], the ratio Pi/r2 is not altered, and the momentum distributions of primary 



nucleons are very similar to the previous calculation. So, the conclusion drawn in ref. [17 
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still holds: a comparision of the calculated spectra with the experimental one favours T n /Y p 
ratios around 2^-3 (or higher), in disagreement with the OPE predictions. On the other 
hand, we have to recall that for a clean experimental extraction of the T n /T p ratio it is very 
important to identify the nucleons which come out from the different non-mesonic processes 
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APPENDIX: SPIN-ISOSPIN NN AND AN -► NN INTERACTIONS 

In this appendix we show how the repulsive A^ A" and AiV Short Range Correlations (SRC) 
are implemented in the A^ A" — > NN and AA" — > A^ A" interactions. The process A^ A" — > A^ A" 
can be described through an effective potentials given by: 

G(r) = g(r)V(r). (Al) 

Here g(r) is a 2-body correlation function, which vanishes as r — > and goes to 1 as r — > oo, 
while V(r) is a meson exchange potential which in our case contains n and p exchange: 
V = V w + V p . A practical form for g(r) is |25[| : 



g{r) = l-j {q e r), (A2) 

With q c ~ 780 MeV one get a good reproduction of realistic A^ A" correlation functions ob- 
tained from G-matrix calculations. The inverse of q c is indicative of the hard core radius of 
the interaction. Since there are not experimental indications, the same correlation momen- 
tum we use for the AA^ interaction. On the other hand, we remind that q c is not necessarily 
the same in the two cases, given the different nature of the repulsive forces. Using the 
correlation function (|A2|) it is easy to get the effective interaction, eq. (|A1|) , in momentum 
space. It reads: 

P 
G NN (q) = K-(g) + V p (q) + -^ {gL(q)qiqj + 9r{q){Sij ~ Mj)} ^i^jT ■ r, (A3) 

where the SRC are embodied in the correlation functions gi and gx- The spin-isospin 
A^Af — > NN interaction can be separated into a spin-longitudinal and a spin-transverse 
parts, as follows: 

G NN (q) = {V L (q)qiqj + V T (q)(S ij - q^)} a^T • r (q t = q t / \ q |), (A4) 



where 



fl 

J IT 



V L (q) = ^{q 2 F 2 n ( q )Gl( q )+g L ( q )} 
mi v > 



(A5) 



fl 



V T (q) = 4 {q 2 C p FZ(q)G p (q) + g T (q)} . 



(A6) 



In the above F p is the pNN form factor (eq. fl2.6| ) with cut-off A p = 2.5 GeV), and G° = 

l/(?o ~ Q 2 ~ m p) i s the P f ree propagator. 

The AN — » A/" A/" effective interactions splits into a p-wave (again longitudinal and trans- 
verse) part: 



G AN ^ NN (q) = {PL(q)qiqj + Pr(q)($ij ~ <Mj)} ^ 



T ■ T, 



with: 



PM = — — {q 2 F*(q)Gl(q) + g£(q)} 



(A7) 



(A8) 



Mq) = ——9t(q), 

m % m % 



and an s-wave part: 



S(q) = ^S{Fl{q)Gl{q)-Fl{ q )Gl{q)} \ q 



(A9) 



(A10) 



Form factors and propagators with a tilde imply that they are calculated changing q 2 — > q 2 + 
7 C - Cp 



g 2 . C p is given by the expression: 



f 2 
r> — j p 

U P ~ 2 

m 



f 2 

J IT 



mi 



The expressions for the correlation functions are the following: 

9L(q) = ~{<1 2 + \<lc} F 2 M)Gl{q) - lq 2 c C p F 2 p (q)G° p (q), 



(All) 



(A12) 



9r(q) = -\qlF 2 M)Gl{q) - [q 2 + \q\ \ C^G^q)., 



(A13) 



9t(q) 



q 2 + $\ F 2 MG Q M\ 



(A14) 
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fl£(?) = -\qIF 2 MG°M. (A15) 

Using the set of parameters: 

q c = 780 MeV, A* = 1.2 GeV, A p = 2.5 GeV, /^/4tt = 0.08, C p = 2, (A16) 

at zero energy and momentum we have: 

9l(0) = g T (0) = 0.615, <?£(0) = <$(0) = 0.155. (A17) 

Howhever we wish to keep the zero energy and momentum limit of ql,t and g^ T as free 
parameters, thus we replace the previous functions by: 

/ x i9l,t{q) a / \ f 9l,t(q) [\-\q\ 

H - T{q) ■" 3 ^m • 3lM 3a ^Wy (A18) 
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